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Abstract. We give integrable quad equations which are multi-quadratic (degree- 
two) counterparts of the well-known multi-affine (degree-one) equations classified by 
Adlcr, Bobcnko and Suris (ABS). These multi-quadratic equations define multi- valued 
evolution from initial data, but our construction is based on the hypothesis that 
discriminants of the defining polynomial factorise in a particular way that allows 
to reformulate the equation as a single-valued system. Such reformulation comes 
at the cost of introducing auxiliary (edge) variables and augmenting the initial 
data. Like the multi-affine equations listed by ABS, these new models are consistent 
in multidimensions. We clarify their relationship with the ABS list by obtaining 
Backlund transformations connecting all but the primary multi-quadratic model back 
to equations from the multi-affine class. 



1. Introduction 

We consider quad equations denned in terms of a polynomial in four variables, Q, that 
is equations of the form 

Q(u, u, u, u) = 0, (1) 

where in the simplest setting u = u(n,m), u = u(n + l,m), u = u(n,m + 1) and 
u = u(n + l,m + 1) are values of a dependent variable taking values in C U {oo} as a 
function of independent variables n, m G Z. The quad equation is called multi-affine if 
the defining polynomial, Q, is degree one in each variable, and we call it multi-quadratic 
if Q is degree two in each variable. 

An important integrability feature that is possible for the quad equation is the 
multidimensional consistency [1, 2]; it has proven to be a natural property of many 
integrable equations in the multi-affine class. Significant work on this is the list of multi- 
affine quad equations with the consistency property obtained by Adler, Bobenko and 
Suris (ABS) in [3, 4]. Quad equations beyond the multi-affine class were considered in [5] 
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where several multidimensionally consistent examples (also beyond the multi-quadratic) 
were obtained. In fact these equations appeared naturally in relation to underlying 
models in a different class, namely the Yang-Baxter maps [6, 7], and emerge from a 
generalisation of results connecting the Yang-Baxter maps with the multi-affine quad 
equations [3, 8, 9]. One multidimensionally consistent multi-quadratic quad equation 
identified in [5] was known earlier due to observations of Adler and Veselov in [10], it 
in fact arises as the superposition principle for Backlund transformations of the KdV 
equation. Recently it has also come to light that the well-known discrete version of 
the KdV equation due to Hirota [11], which is multi-affine but absent from the ABS 
list, is naturally understood within the consistency framework as a special case of a 
multi-quadratic quad equation [12]. The few known examples therefore indicate that 
multi-quadratic quad equations are a quite natural and potentially rich class of integrable 
systems. 

A feature of higher degree discrete models, like the multi-quadratic quad equations, 
is that they define multivalued evolution from initial data. This adds richness, but 
also another level of difficulty in dealing with such systems. This difficulty is however 
mitigated in the mentioned integrable cases because they can be reformulated as a single- 
valued system with augmented initial data. For the models in [5] this is because the 
variables present in the associated Yang-Baxter map themselves satisfy a single-valued 
system. But that situation can actually be viewed as a special case of a more general 
constructive approach to this kind of reformulation [12]. In this approach auxiliary 
variables are introduced on lattice edges, which are similar to variables of an associated 
Yang-Baxter map, however, rather than satisfying an independent system, they instead 
participate in a mixed system involving both vertex and edge variables. The resulting 
model is of a type similar to the class introduced in [13] but with the additional feature 
of preserving algebraic relations on the lattice edges. The reformulation procedure relies 
completely on a discriminant factorisation property of the defining polynomial. It is this 
special property that is the departure point and the main emphasis of the present work, 
we use it as a foundation to enable some systematic investigations of the multi-quadratic 
quad equations. 

The main result presented here is a list of multi-quadratic quad equations with the 
discriminant factorisation property. To construct the models we start from biquadratic 
polynomials (the factors of the discriminant) which are associated with the edges of the 
lattice, and this provides a natural correspondence with the multi-affine ABS equations 
whose construction also involves edge biquadratics. Due to the factorised discriminant 
hypothesis the models we list admit reformulation as single-valued systems. However 
the more important property of these equations, which is actually rather remarkable 
because it is not explicitly built into the construction, is the integrability feature of 
multidimensional consistency. 

The second part of our paper is devoted to developing the transformation theory 
of the obtained models. The most pressing question we seek to address relates to 
existence of transformations of Backlund or Miura type connecting these multi-quadratic 
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models back to the better studied multi-affine ABS equations. Such transformations 
are not inherent from our construction, rather we tackle the problem of obtaining 
transformations a posteriori using separate methods. The main result is to connect 
all except the primary model, namely the multi-quadratic counterpart of Q4, back to 
equations in the multi-affine (ABS) class. 

We proceed as follows. In Section 2 we explain the basic discriminant factorisation 
hypothesis for the multi-quadratic quad equations. Additional assumptions are 
explained in Section 3 and the obtained models are listed in Section 4. The method to 
reformulate these models as systems that define single- valued evolution from initial data 
is given in Section 5 and their multidimensional consistency is described in Section 6. 
Section 7 is devoted to the transformation theory of the models, in particular we recall 
the multi-affine equations from the ABS list in Section 7.1 and give transformations 
connecting the multi-quadratic models to them in Section 7.2. The methods used to 
obtain these transformations are explained in Sections 7.3, 7.4 and 7.5. Some questions 
raised by the results reported here are discussed in Section 8. 

2. The factorised-discriminant hypothesis 

In the multi-quadratic case considered here, the defining polynomial Q of the quad 
equation (1) is of degree two in each of the four variables. To calculate the evolution 
defined by (1) from some initial data requires solving this equation locally as a quadratic 
equation for one of the arguments, and in particular taking the square root of its 
discriminant with respect to that argument. In general this discriminant is a polynomial 
of degree four in each of the remaining three arguments; the property of polynomial Q 
that we study here is the factorisation of this discriminant into a product of three factors: 



where Hi and H 2 are degree-two polynomials in each of two variables, i.e., biquadratics, 
and Gi, . . . , G 4 are degenerate biquadratics of the form Gi(u, u) = (ai + biU+CiU+diUu) 2 , 
i G {1, 2, 3, 4}, i.e., the square of a polynomial which is degree-one in each variable. This 
is the most general hypothesis that allows replacement of (1) by an associated single- 
valued system through introduction of auxiliary edge variables <J\ and a 2 satisfying the 
edge relations 



The key features allowing this are that all discriminants in (2) are squares which clearly 
leads locally to a rational model, and furthermore that the edge relations on opposite 
sides of the quad are the same, allowing the rational model to replace the quad equation 
globally (for instance in the simplest setting throughout Z 2 ). The procedure to obtain 



A[Q(u,u,u,u),u] 
A[Q(u,u,u,u),u] 
A[Q(u,u,u,u),u] 
A[Q(u,u,u,u),u] 



Hi (u, u)H 2 (u, u)Gi (u, it) , 
Hi(u, u)H 2 (u, u)G 2 (u, u), 
Hi(u, u)H 2 (u, u)G 3 (u, u), 
Hi(u, u)H 2 (u, -u)G 4 (-u, u), 



(2) 



o\ = Hi(u, u), o\ = H 2 (u, u). 



(3) 
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the single-valued system is straightforward, it was explained in detail in [12] and an 
example will also be included later in this paper (Section 5). 

We remark that the restriction to Gi, . . . , G4 being degenerate seems unnatural in 
some regards, in fact there exist sets of polynomials {Q, Hi, H 2 , H 3 } satisfying system 
(2) with Gi = G 2 = G 3 = G 4 = H 3 , and where H 3 is also non-degenerate. These 
more symmetric solutions of the problem are interesting, however to understand the 
sense in which the resulting polynomial Q defines an integrable discrete model requires 
a substantial alteration of the setting, and this will be studied in detail elsewhere. 
Also we remark that the assumed symmetry can be relaxed, for instance replacing Hi 
appearing in the third and fourth equations of (2) with Hi ^ Hi. Such systems are of 
interest too but the associated polynomial Q in that case is more naturally interpreted 
as defining a Backlund transformation between models. 

3. Additional assumptions 

The solution of (2), in the sense of obtaining a set of polynomials 



for which system (2) is identically satisfied, is a difficult classification problem in its full 
generality. Further assumptions allow this problem to be solved using computer algebra, 
these come from looking at examples identified in [12]. 

We assume invariance of Q when the arguments are permuted by the Klein four- 
group (the Kleinian symmetry) 



{Q, Hi, H 2 , Gi, G 2 , G 3 , G 4 } 



Q(u,u,u,u) = Q(u,u,u,u) = Q(u,u,u,u). 



(4) 



This reduces the generic multi-quadratic polynomial Q to the form 




(5) 
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for some set of coefficients {ai, . . . , 027}. This symmetry is a strong additional 
assumption, for instance by itself it is actually sufficient for integrability of quad 
equations in the multi-affine class [4, 14]. 

We also assume that both biquadratic polynomials Hi and H2 are symmetric, and 
furthermore that they are taken from one of the following one-parameter families (the 
parameter being denoted by p) 

q4* : ^ (l + u 2 u 2 + u 2 p 2 + p 2 u 2 ) - (u 2 u 2 p 2 + u 2 + u 2 + p 2 ) - ( c 2 - \ j uup, (6) 

(7) 

(8) 

(9) 

(10) 

(11) 
(12) 

(13) 

Where it appears 5 G {0, 1}, whilst in (6) c G C \ {0, 1, — i} is a fixed constant. 
Up to a point transformation of the parameters these biquadratic polynomials coincide 
with those associated with the ABS equations [3, 4] (precise parameter associations will 
be given later in Section 7.1). Note that (12) and (13) are considered to be biquadratics 
here also, but this is in a projective sense, so loosely speaking we view for instance (13) 
as the polynomial {u + u + p) {u — 00) {u — 00). We remark that by a Mobius change of 
variables a symmetric biquadratic polynomial can always be brought to one of the forms 
(6)-(13) for some choice of parameter p, or else is the square of a multi-affine polynomial 
(like Gi, . . . , G4 above), a possibility which we therefore exclude by restricting to (6)- 
(13). The discriminant of biquadratics (6)-(13) with respect to u is a polynomial of 
degree at most four in u, the main characterising feature of the biquadratic families is 
that roots of this discriminant polynomial do not change upon altering parameter p. 
Finally, the assumed symmetry of Q together with (2) implies already that 

1 — (_J2 — <-^3 — ( - r 4 — • 

in fact we will go further and make the assumption that 

G(u,u) = {u-u) 2 , (14) 
which is our last additional assumption. 



q3* : 




1) — uup, 


q2* : 


1 ^ + u 2 +p 2 )- \ 


{uu + up + pu) , 


ql* : 


[u-uf -p, 




a2* : 


-(1 + u 2 u 2 ) — puu, 
2 




al* : 


(u + u) 2 - p, 




h3* : 


uup + 5 2 , 




h2* : 


u + u + p. 
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4. The obtained multi- quadratic quad equations 

The additional assumptions enable solution of the discriminant factorisation hypothesis 
directly by computer algebra, which gives the main result of our paper as follows. 

Proposition 1 Let Q be a polynomial of degree two in each of four variables with the 
Kleinian symmetry 

Q(u,u,u,u) = Q(u,u,u,u) = Q(u,u,u,u), (15) 

such that the discriminants factorise as follows 

A[Q(u, u, u, u) , u] oc (u — u) 2 Hi(u, u)H 2 (u, it), (16) 

where Hi and H 2 are biquadratic polynomials. If Hi and H 2 are taken from one 
of the biquadratic families (6)-(13) with generic parameter choices denoted p and q 
respectively, then Q is determined up to an overall constant and defines respectively the 
quad equations (17)-(24) below. 

Q4*: 

(p — q)[(c~ 2 p — c 2 q)(uu — uu) 2 — (c~ 2 q — c 2 p)(uu — uu) 2 } 

-(p - q) 2 [(u + u) 2 {\ + u 2 u 2 ) + {u + u) 2 {\ + u 2 ti 2 )} ^ 
+[(u — u)(u — m)(c~ 1 — cpq) — 2{p — q)(l + uuuu)} 
x [{u — u)(u — u)(c -1 p5 — c) — 2(p — q){uu + uu)] = 0, 

Q3*: 

(p — q)[p(uu — uu) 2 — q(uu — uu) 2 ] — 5 2 (p — q) 2 [{u + u) 2 + (u + it) 2 ] 
+[(u — u){u — u) —25 (p — q)][(u — u){u — u){pq — 1) — 2{p — q){uu + uu)] = 0, 
Q2*: 

(p — q) [p(uu — uu) (u + u — u — u) — q{uu — uu) (u + u — u — u)] 

+(u — u)(u — u)\p(u — u){u — u) — q{u — u)(u — u) — pq{p — q)] = 0, 

Ql*: 

(p - q)[p(u + u - u - u) 2 -q(u + u-u- u) 2 } 

+A(u — u)(u — u)[p{u — u)(u — u) — q(u — u)(u — u)] = 0, 

A2*: 

(p - q) [p(uu - uu) 2 - q(uu - mm) 2 ] 

+(u — u)(u — u)[(u — u)(u — u){pq — 1) + 2{p — q)(l + uuuu)} = 0, 

Al*: 

(p — q) [p{u — u + u — u) 2 — q(u — u + u — u) 2 ] 

+4(-u — it) (u — it) [p{u + u)(u + u) — q{u + u)(u + u)] = 0, 

H3*: 

(p - q) [p(uu - uu) 2 - q(uu - uu) 2 } 

+(u — u)(u — u)[(u — u)(u — u)pq — 45 2 (p — q)] = 0, 
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H2*: 

(P - Q)[p( u -u + u-u) 2 -q(u-u + u- uf] 

+(u — u)(u — u)[(u — u){u — u) — 2{p — q)(u + u + u + u)] — 0. 

We have denoted the list of models here QA*,...,H2* due to their natural 
correspondence with equations from the ABS list, which are denoted QA, . . . ,H1 [3]. 
(The ABS list and its precise relationship with (17)-(24) will be given in Section 7.1.) 
The biquadratic polynomials associated with ABS equations Qls=o, ^41<5=o and HI are 
already squares, which means the multi-quadratic counterparts of those models factorise 
into products of multi-affine polynomials, cases we exclude here as degenerate. 

Similar to the situation for ABS equations, all of (18)-(24) can be obtained 
from the primary model (17) by limiting procedures. All listed models are non- 
equivalent by autonomous Mobius changes of variables and point transformations of the 
parameters. However, the models A2* and Q3* s=0 are related by the non-autonomous 
point transformation u — > •u( _1 ) n+m , whilst Al* and Ql* are related by u — >■ (— l) n+m u. 

To our knowledge the equations in the list (17)-(24) are new except for the following 
cases. The model H2* (24) appeared as the superposition principle for solutions of the 
KdV equation in [10] and in relation to Yang-Baxter maps in [5]. The model Al* (22) 
was obtained, although not in explicit form, in [15]. The model A2* (21) was obtained 
originally in [12] as the superposition principle for Backlund transformations of Hirota's 
discrete KdV equation [11]. 



5. Reformulation as single- valued systems 

One of the most salient features of quad equations from the multi-quadratic class 
is that they define multivalued evolution from initial data. However, built into the 
construction of models (17)-(24) is the discriminant factorisation property, and this 
allows the quad equation to be reformulated as a system that defines single-valued 
evolution from augmented initial data [12]. Here we perform this reformulation on the 
primary model Q4* (17), and describe the relationship between the quad equation and 
its reformulation. 

As described in Section 2, the auxiliary variables on lattice edges (see Figure 1) are 
introduced through the relations 



uup, 

(25) 

uuq, 



o* = ^(l + u 2 u 2 + u 2 p 2 + p 2 u 2 ) - ^- (u 2 u 2 p 2 + u 2 + u 2 + p 2 ) - (c 2 - ^\ 

a\ = ^ (l + u 2 u 2 + u 2 q 2 + q 2 u 2 ) - ^- [u 2 u 2 q 2 + u 2 + u 2 + q 2 ) - ^c 2 - ^ 

which are in terms of the associated biquadratic polynomial given earlier in (6). By 
solving equation (17) as a quadratic equation for u and exploiting these edge variables 
we obtain an equation of the form 

J-(u, u, u, u, <J\<J2) = (26) 
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u. 



o-i 



,u 



0-2 



0-2 



u 



0"! 



u 



Figure 1. Variables assigned to the vertices of a quadrilateral, and auxiliary variables 
to the edges. 



(28) 



which is by construction of polynomial degree one in u and o"i<r 2 , for instance: 
J^iu, u, ii, u, 0\02) := u\[cT 2 p — c 2 q)(uu — uu) + (c~ 2 q — c 2 p)(uu — uu)} 

— {u — -u)[(c _1 — cpq)(u 2 + uu) + (c~ 1 pq — c)(l + u 2 uu) + 2o" 1 o" 2 ] (27) 

— (p — q)(u — u)(l + u 3 u) 

(the precise form of this expression is chosen because it is simple, but it is not unique). 
Sequentially solving for each of the quad variables of (17) in the same way leads to the 
following system 

^(u, -u, u, u, aia 2 ) = 0, 
F(u, it, m, u, o-icr 2 ) = 0, 

J-"(-U, U, U, U, 0"iO" 2 ) = 0, 

J-"(-u, u, u, u, <5"i5" 2 ) = 0, 

up to the choice of sign of the discriminant terms appearing in the last argument. 

The system (28), (25) is the single-valued model associated with (17). It is easily 
verifiable that each equation in (28) implies (17) is satisfied modulo the relations 
(25). The signs of the discriminant terms in (28) have been chosen for self-consistency 
(ensuring that any one equation in (28) is a consequence of the other three), to preserve 
the Kleinian symmetry, and to preserve the consistency property described in the 
following section. 

The usual initial value problem for a quad-equation involves specifying the 
dependent variable on vertices along some admissible lattice path or collection of paths. 
The modification required for system (28) is that variables on path edges, subject to the 
constraint (25), also need to be specified. The system (28), (25) is therefore actually 
a model of vertex-bond type as described in [13], but it has the additional feature of 
preserving algebraic relations (25) on lattice edges. In this setting the multivalued 
evolution defined by equation (17) is reflected in the non-uniqueness of initial edge 
variables when they are defined in terms of the initial vertex variables through (25). 

In Table 1 the polynomial J 7 that appears in the single-valued system associated 
with each equation from the list (17)-(24) is also given. The auxiliary variables present 



Multi- quadratic quad equations: integrable cases 



9 





M[(C p - 


- c q)[uu — uu) + (c q — c p)[uu — uu)\ 




-(«- 


u)[(c~ l — cpq){u 2 + tra) + (c~ l pq — c)(l + u 2 uu) + 2<7i<7 2 ] 




-(p- 


g)(« - m)(1 + w 3 -u) 


no* 


m[p(mm — 


mm) + g(-uw — mm)] — (m — m)[m 2 + mm — 5 2 (1 — pg) — 2(7! er 2 ] 






-g)(w-w) 


Liz 


m[p(m + m — m — m) + — m + m — u)] + p{uu — mm) + q(uu — uu) 




-(«" 


m) [(m — u) (u — U)+pq + 4<7i<72] 


0,1 


p(m + m - 


- M — M) + g(M — M + M — M) — 2(m — M)[(M — m)(m — m) — 0"iO"2] 


A2* 


m[j>(mm — 


mm) + q{uu — mm)] — (m — m)(1 + m 2 mm — 2<7i(7 2 ) 


Al* 


p(m + m - 


- M — m) + ?(m — M + U — u) — 2(m — m)[(m + m)(m + U) — 0"iO" 2 ] 


H3* 


u[p{uu — 


mm) + q(uu — uu)] + 2(m — u)(5 — O\0~-2) 


H2* 


p(u + u - 


- u — u) + q[u — u + m — u) + (m — u) (m + m + 2m — 2<t 1 o- 2 ) 



Table 1. Polynomials u, u, m, 0102) that define the single-valued system (28) 
associated with the multi-quadratic models (17)-(24) of Proposition 1. The auxiliary 
variables u\ and cr 2 satisfy edge relations a\ = H\{u, u) and o\ = H 2 (u, u) in terms of 
the associated biquadratic polynomials listed in (6)-(13). 



in the table are introduced through the edge relations which can be written generically 

as 

(t\ = H 1 {u,u), a 2 2 = H 2 (u,u), (29) 

in terms of the associated biquadratic polynomials (cf. Proposition 1). 

The single- valued reformulation (28), (29) that we have proposed for the quad 
equation (1) is autonomous, but there is the possibility to make a non- autonomous 
reformulation of the model in which signs of discriminant terms in (28) change from quad 
to quad. Such non-autonomous reformulation is more difficult to study, however in the 
autonomous case the quad equation and its reformulation are not globally equivalent. 
They are of course locally equivalent, by which we mean that on a single quad the 
equation (1) is sufficient for consistency of system (28) in o~i, <r 2 , allowing the edge 
variables to be constructed from vertex variables on a single quad. The simplest setting 
in which the quad equation and its autonomous reformulation are not equivalent is a 
square domain involving nine points (four quads). On this larger domain we can write 
the system 

.F(m, m, m, m, (X1O2) = 0, 

J r (u,u,u,u,g 1 a 2 ) = 0, ^ 
J-(u : M, M, M, o~iq 2 ) = 0, 
J-(u, M, M, M, 0" 1 cr 2 ) = 0. 

where u = u(n — l, m), u — u(n,m — l) etc, this system being a consequence of imposing 
(28). The equation (1) imposed on each quad is not sufficient for consistency of (30), 
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specifically elimination of o\,ai,q x and q 2 from (30) yields a constraint on the nine 
vertex- variables which is not a consequence of (1). In particular (30) can be used to 
obtain u uniquely from data {u, u, u, u, u, u, u, u}, whereas if (1) alone were imposed this 
data would yield two possible values of u. This example therefore demonstrates that 
global equivalence of the quad equation and its associated single-valued system is not 
possible unless a non-autonomous reformulation would be considered. 

More generally than the example above, if initial data is specified on verticies along 
some admissable lattice path, then the autonomous single- valued system (28), (29) 
generically determines 2 < ^ 1 ) solutions where x is the number of edges along the path 
(corresponding to the number of auxiliary edge variables), whereas using the equation 
(1) alone the number of solutions determined would be 2 y where y is the number of 
quads in the domain. Because (on a finite domain) y > x — 1, the autonomous single- 
valued reformulation of the system can be said to limit the multivaluedness from initial 
data. 

6. Multidimensional consistency 

The purpose of this section is to explain a key integrability feature of models (17)- 
(24), namely their multidimensional consistency [1, 2, 3]. This emerges quite naturally, 
although it has not been explicitly built into the construction of these equations. 

The multidimensional consistency involves not just one equation in isolation, and 
in fact for the equations listed it involves the whole family of equations obtained by 
varying parameters p and q. The consistency is between members of this family with 
different choices of the parameters. Due to the central role they play it is therefore 
convenient to explicitly include the dependence on parameters p and q when writing a 
generic equation from (17)-(24): 

Q P; q(u,u,u,u) = 0. (31) 

The key properties (involving the parameter dependence) of the generic defining 
polynomial are first the symmetry 

Q p , q (u, u, it, u) = Q q>p (u, it, u, u), (32) 

and second the consistency of the system 

Qp,q(u, U, U, u) = 0, Qp, g (u, U, U, u) = 0, 

Q q ,r{u, it, it, it) = 0, Qq, r (u, U, U, it) = 0, (33) 
Qr,p(u, U, U, U) = 0, Q r ,p(u, U, U, U) = 0, 

for any choice of the parameters p, q and r, which is usually visualised by assigning 
variables to vertices of a cube as in Figure 2, and equations to faces. By consistency we 
mean that for generic initial data {u, u, u, u} the system (33) has at least one solution 
for the remaining variables {u,u,u,u\. Directly by polynomial manipulation it can be 
confirmed that for each quad-equation (17)-(24) the system (33) is consistent in this 
sense, and in fact determines four possible solutions {u, u, u, u}. 
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Figure 2. Variables assigned to the vertices of a cube. In the case of the multi- 
quadratic quad equations (17)-(24) the system (33) determines four possible values of 
u from initial data u, u, u and u. 



Alternatively, and more explicitly, the consistency property can be verified by 
reformulating the equation as a single- valued system as described in Section 5. For 
instance in the case of the primary model (17) a direct calculation yields the following 
expression that determines u in terms of the initial data on the cube: 

u(c 2 — c~ 2 )[p(q — r) 2 (uu — uu) + q(r — p) 2 (uu — uu) + r{p — q) 2 {uu — uu)] 

— {q — r)(r — p){u — -u)[(c _1 — cpq){u 2 + uu) + (c~ x pq — c)(l + u 2 uu) + 2<7i<7 2 ] 

~ _ _ ( 34 ) 

— (r — p)(p — q){u — u)[(c 1 — cqr){u 2 + uu) + (c 1 qr — c)(l + u 2 uu) + 2g20~z\ 

— {p — q)(q — r)(u — ■u)[(c~ 1 — crp){u 2 + uu) + (c _1 rp — c)(l + u 2 uu) + 2<t 3 0"i] = 0, 
where a±, a 2 and gt 3 are determined from the initial data through equations (25) and 

o-l = - (1 + u 2 u 2 + u 2 r 2 + r 2 u 2 ) - — (u 2 u 2 r 2 + u 2 + u 2 + r 2 ) - ( c 2 - \ j uur. (35) 

Although the transformation group (<7i, o"2, o" 3 ) i— )■ (±o"i, ±o" 2 , ±^"3) leaving (25), (35) 
invariant contains eight elements, the initial data {u, u, u, u} leads to only four distinct 
values of u due to the symmetry (<Ji, o- 2 , <t 3 ) i-)- (— a±, — o~2, — cr 3 ) of (34). 

It is important to note that the consistency property of the associated single- valued 
system can be broken by reversing sign of all discriminant terms in (28). This subtlety 
reflects the fact that polynomial system (33) determines only four possible values of u 
from the initial data, and not eight as might be expected. Other polynomials given in 
Table 1 are also chosen with the consistency property in mind, so that in particular the 
sign chosen for the discriminant terms is important. 

The consistency property of models (17)-(24) means that their more general setting 
is for a dependent variable u defined on Z d , where d > 1 is the dimension, determined 
by the system 

Q Pi , Pj (u, JiU, Tju, JJju) = 0, i, j e {1 . . . d}. (36) 

Here pi, . . . , pd are some set of parameters, while T 1; . . . , T^ are shift operators in each 
dimension. Probably the simplest Cauchy problem for (36) is to specify initial data on 
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coordinate axes, but actually this multidimensional setting is the departure point for 
an extremely rich variety of initial value problems and lattice configurations [2, 10, 16]. 
This also aligns with the notion of solvability for models in statistical mechanics [17]. 
From a slightly different perspective the consistency property yields a great deal of 
control over the solution structure of these models by providing immediately a natural 
auto-Backlund transformation. This allows for example to construct exact solutions 
as periodic Backlund chains (the discrete analogue of finite-gap solutions), as well as 
soliton-type solutions from Backlund iteration. 

We remark that on the domain of a hypercube, {0, l} d , associated edge variables 
can be obtained (as described in Section 5) for any solution of (36), this domain is 
therefore singled out as one on which the quad equation and its autonomous single- 
valued reformulation are globally equivalent. 

7. Transformations to ABS equations 

All equations listed by ABS in [3], with the exception of Q4, are mutually related 
via Backlund or Miura type transformations [18, 19, 20, 21]. Table 2 here lists 
similar transformations connecting all except the primary multi-quadratic model Q4*, 
specifically equations (18)-(24), back to multi-afline equations from the ABS list. Before 
describing these transformations in more detail we first recall the ABS list in (37)-(45) 
below. 



7.1. The ABS list 



Q4: 



sn(a)(iw + vv) — sn( ( 5)(v'0 + vv) 

— sn(o; — (3) [vv + v v — k sn(a) sn(/3) (1 + vvvv)] = 0, 



(37) 



Q3: 



(a 



l/a)(vv + vv) — ((3 — 1/ j3){vv + vv) 

-(a//3 - f3/a)[vv + vv + 5 2 (a - l/a)(/3 - l//3)/4] = 0, 



(38) 



Q2: 



a(v 



v)(y — v ) — (3{v — v)(v — v) 

+a/3(a - fi)(v + v + v + v - a 2 + a/3 - f3 2 ) = 0, 



(39) 



Ql: 



a(v — v)(v — v ) — (3{v — v)(v — v) + 5 2 a(3(a — (3) = 0, 



(40) 



A2: 



(a 



l/a)(vv + vv) -{(3- l/(3)(vv + vv) - (a/ (3 - (3/a)(l + vvvv) = 0, 



(41) 



Al: 



a(v + v)(v + v) - (3(v + v)(v + v) - 5 2 a(3(a - (3) = 0, 



(42) 
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H3: 



H2: 



HI: 



a(vv + vv) — p{vv + vv) + 5(a 2 — P 2 ) = 0, 



(v — v)(v — v) — (a — j3)(v + v + v + v + a + P) — 0, 



(43) 



(44) 



(v-v)(v-v)-a + p = 0. (45) 

In (37)-(45) we have reproduced the ABS list. Where it appears 5 G {0, 1} and 
for Q4 (37) k e C \ {0, 1, —1} is the modulus of the Jacobi elliptic function sn. The 
parametrisations coincide with the ones given by ABS in [3] except for the case of Q4, 
this canonical form (the Jacobi form) was obtained by Hietarinta in [22]. Each ABS 
equation is defined by a polynomial Q = Q(v, v, v, v) that is degree one in four variables. 
It was shown in [3] that this polynomial is characterised in terms of biquadratics through 
a generalised discriminant formula as follows 

(dcQW-Q) - (d d d~Q)Q kH^v), (d v Q)(d,Q) - (d v d,Q)Q oc H^vj), 

{dvQ)(d~Q) - (dvd~Q)Q oc H 2 (v,v), (d v Q)(d,Q) - (d v d,Q)Q oc H 2 {v,v). 

The polynomials Hi and H 2 appearing in (46) coincide with those in (16) up to an 
overall constant, provided we make the following associations between parameters a and 
P appearing in the multi-afnne equations (37)-(44) and parameters p and q appearing 
in their multi-quadratic counterparts (17)-(24), 



(46) 



Q4 
Q3,A2 

Q2 
Q1,A1 

H3 

H2 



y/k sn(a + K), 
(a + l/a)/2, 



q = Vksn{p + K), 
q=(P + l/P)/2, 

q = P\ 
q = 5 2 P\ 
q = l/P, 
q = 0, 



Vk, 



(47) 



p 
p 

p = a-, 
p = 5 2 a 2 , 

p = 1/a, q = l/P, 6 2 -+ 5, 

p = a, 

where in the case of Q4 the parameter K is standard notation for the quarter period of 
the sn function (with modulus k) satisfying sn(i^) = 1, sn'(K) = 0. 

We remark that for the multi- afline equations HI (the lattice potential KdV 
equation [23, 18]), Ql<5=o (the lattice Schwarzian KdV equation [18]) and Al5 =0 the 
associated biquadratic is the square of a polynomial which is degree one in each variable, 
and as mentioned before their multi-quadratic counterparts factorise into products of 
multi-affine polynomials which we have excluded as degenerate cases here. 



7.2. An example transformation 

To explain precisely the meaning of the entries in Table 2 we give here an explicit 
example, specifically the second entry of the table. Consider the coupled system of 
equations 

(v — v)(uv — uv) + avv = 0, (v — v)(uv — uv) + Pvv = 0, (48) 
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Eq. in u Backlund transformation Eq. in v 

Q3; =2a2 _ 1 a[(u + 5)v 2 + {u + 5)v 2 } = (u + u + 25a 2 )vv Q3 5=0 

Q2* =Q (v — v)(uv — uv) + avv — Ql«s=o 

A2 ;= 2 ( Q2+ i)/(a^i) ®(v 2 v 2 + 1) = vv[uu(a 2 - 1) + (a 2 + 1)] A2 

H3* =4/q2 uu = vv(vv + 5 a) H3 

H3; =4/Q 4auu = ({; + v) 2 - 5 2 a 2 Al 

A1* =q2 2vv(u + u) = a{v 2 v 2 + 1) H3 



<5=0 



Ai; =Q 2(v + v) (u + u) = (v + v) 2 + a Al 



<5=0 



H2* _ a u + u-a = {v + v) 2 HI 



A2 p=2a 2 -1 


a(uuv 2 v 2 + 1) = (uu + l)vv 


A2 


H3* =4/a 


Aauu — (v — v) 2 — 8 2 a 2 


Ql 


Ql* 2 


2vv(u — u) — a(v 2 v 2 + 1) 


H3^=o 


qi;= q 


2(v — v)(u — u) = (v — v) 2 + a 


Q 1-5=0 



Table 2. Backlund transformations connecting equations from the multi-quadratic 
class to the multi-affinc (ABS) class. Details of implementing the transformations 
are explained in Section 7.2. Transformations are given up to composition with point 
symmetries of the equations in u and v. For completeness we include in the second part 
of the table some transformations that can be obtained from those in the first part 
by composition with (non-autonomous) point transformations. The transformation 
connecting H2* to HI was given originally in [10] up to composition with point 
symmetries. Methods to obtain the transformations are described in the text, first 
by non-symmetric degeneration of the auto-Backlund transformation in Section 7.3, 
second by taking advantage of a natural connection with Yang-Baxter maps in Section 
7.4, and third in Section 7.5 by exploiting the fact that the transformations, like the 
equations themselves, can be characterised by discriminant properties of the defining 
polynomial. 



which involve two functions u = u(n,m) and v = v(n,m), where as usual n, m e Z, 
u = u(n + 1, m) and u = u(n, m + 1), etc. 

For a fixed function v = v(n,m) the equations (48) are coupled discrete Riccati 
equations for u, which are compatible if v satisfies Ql<5=o. Due this choice of v the 
function u = u(n,m) that emerges as the solution of (48) then satisfies equation Q2*, 
that is (19), with parameter associations p = a and q = (3 (notice that this parameter 
association is different from the one in (47)). For brevity the second equation in (48) 
and the parameter association between q and (3 are omitted from Table 2 because they 
can be inferred from the first equation in (48) and the association between p and a that 
have been listed. 

The transformation defined by (48) from v to u is therefore a quite standard 
discrete- Riccati type of Backlund transformation, see for instance [18, 2, 19]. However 
the system (48) also defines an inverse transformation from uto v which is less standard 
because the equations for v are not of Riccati type. Similar to verifying the consistency 
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property of the multi-quadratic models (17)-(24) (cf. Section 6), it is convenient to 
handle this system by introducing auxiliary variables. Specifically for the example (48) 
the auxiliary variables enter through the edge relations 

er? = - (u 2 + u 2 + a 2 ) (uu + ua + au) , 

i I (49) 

a i = t (« 2 + + P 2 ) - o + u/3 + /3u) , 

so in fact they are the Q2* auxiliary variables. They allow (48) to be re-written as 

2uv = (u + u — a — 2o~i)v, 2uv = (u + u — (3 — 2<72)i>, (50) 

which have degree-one polynomial dependence in v, v and v. Implementation of the 
transformation from u to v therefore requires a solution (u,a 1 ,a 2 ) of the single-valued 
reformulation of Q2* (cf. Table 1). From such solution the system (50) determines the 
function v satisfying Q 1,5=0. 

Similar to the example we have focused on in this section, all transformations listed 
in Table 2 are standard Riccati-type systems for u. The inverse transformation to 
obtain v is always more involved, reducing to a Riccati type system for v only after the 
introduction of auxiliary variables associated with the solution of the multi-quadratic 
quad equation for u. 

7.3. Obtaining the transformations 

It is rare to have a systematic method to obtain non-local transformations of Backlund 
or Miura type between a given pair of equations. The transformations obtained 
here do however fit into a general framework, specifically they are defined in terms 
of polynomials which are characterised by their discriminants, a framework which is 
therefore consistent with the main theme of this paper. We exploit this general point of 
view in Section 7.5 however, we instead use more direct methods to obtain most of the 
transformations listed in Table 2. A method involving non-symmetric degeneration of 
the auto-Backlund transformation is explained in this section, transformations obtained 
through a connection with Yang-Baxter maps will be explained in detail in the following 
section (Section 7.4). 

There is a method to obtain Backlund transformations between distinct quad 
equations that was developed in [19], it is constructive insofar as it takes as 
a starting point the already obtained equations. It relies on the natural auto- 
Backlund transformation, which for the equations in question is inherent from the 
multidimensional consistency, and exploits this in conjunction with the hierarchical 
relationships between equations. Specifically we seek to connect an equation to its 
degenerate counterpart by making a non-symmetric degeneration of the auto-Backlund 
transformation. 

The principal example where the non-symmetric degeneration technique has been 
used is to obtain the second transformation in Table 2 (which was also the example 
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considered in Section 7.2). We exploit the fact that the substitution 

u= , - . (51) 
e(l + v) V ; 

into equation Q2* yields equation Ql«5=o to first order as e — > 0. The natural auto- 
Backlund transformation for Q2* (19) is as follows 

Qp,r(u, U, V, V) = 0, Qq,r(u, U, V, v) = 0, (52) 

where <2 M is the defining polynomial of this equation, and the transformation connects a 
solution u — u(n, m) to another solution of the same equation v — v(n,m). Applying the 
degeneration procedure to solution v and judiciously choosing the Backlund parameter 
we are led to write 

Qp, r (u, u, rv/(l + v),rv/(l + v)) = 0, Q q , r {u, u, rv/(l + v),rv/(l + v)) = 0, (53) 

which at leading order as r = 1/e — > oo yields exactly the desired non-auto Backlund 
transformation (48). 

7.4- From Yang-Baxter maps to Backlund transformations 

The Yang-Baxter maps given in [24], when suitably interpreted as a system of equations 
for functions defined on edges of the lattice, are naturally connected with multi-affme 
quad equations from the ABS list [3], in particular a potential function for the edge 
variables is governed by a quad equation. A systematic method to obtain the Yang- 
Baxter system on edge variables starting from the multi-affme equations from the 
ABS list was given in [9]. Developments reported in [5, 15] go further, but in the 
opposite direction, showing that also non-multi-affme multidimensionally consistent 
quad equations can emerge as potential for the Yang-Baxter systems. And furthermore, 
that two different quad equations emerging in this way from the same Yang-Baxter 
system, as we shall see in this section, can immediately be connected through a 
Backlund-type transformation. 

Of particular relevance here are the models that have been considered in [25] 
alongside the present work. This section is devoted to recalling the relevant models 
from [24] and giving details of the theory developed in [5, 15, 25] which, combined with 
the methods developed in this paper, can be used to obtain many of the transformations 
in Table 2. 

The starting point is a system of equations for two variables, say s and t, assigned 
to the edges of Z 2 oriented in the n and m directions respectively (similar to a\ and 
<t 2 introduced in Section 5, cf. Figure 1). The particular systems relevant here are as 
follows. 
Ff. 

a{l - (3 2 )s - (3{1 - a 2 )t - a 2 + (3 2 

s — t 



0(1 - a 2 )s - a(l - (3 2 )t + (a 2 - (3 2 )sf 

^l^ a 2)t- a (i-(5 2 ) s -f3 2 + a 2 { } 

"ail - f3 2 )t - /3(1 - a 2 )s + (f3 2 - a 2 )sV 
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1 1 '' 



in- 



t 



t as - pt - 5(a 2 - P 2 ) 
ps-at 

as-pt- 5(a 2 - P 2 ) { } 

Ps — at 

as — Pt 

Ps — at ' 

a — P 
s = t+ 

- S a-P < 57 > 

which are nothing but the quadr national Yang-Baxter maps presented in [24] suitably 
interpreted as an equation on the lattice (we have omitted model F IV because it is 
not used here). The relevant feature of systems (54)-(57) is that they admit a three- 
parameter family of potentials, we denote the potential here by /. The potentials 
corresponding to these models, which are derived in [25], are as follows: 

f + f = A ln(s) + B ln(s - a) + C Has - 1) + \{B + C) ln(« 2 - 1), 
/ + / = A\n(t) + B\n(t - 0) + C\n(pt - 1) + \{B + C) ln(/3 2 - 1), 

f + f = Aa(2s-5a) + Bki(s) + CHs-5a), 
[rn) (59) 
/ + / = Ap(2t-Sp) + Bln(t) + Cln(t-6p), 

,j? \ f + f = AHs) + aBs + C/s, 

{Fin) ~ (60) 
f + f = A\n(t) + PBt + C/t, 

f + f = As + B(s 2 + a) + C(s 3 + 3as), 
(Fv) - (61) 
/ + / = At + B(t 2 + p) + C(t 3 + 3pt), 

The parameters A, B and C may be chosen freely in each case. 

For the purpose of illustration we focus on the system (54). In particular writing 
/ = ln(f) for the potential corresponding to the choice of the parameters A — 1, 
B = C = in (58) we find 

vv = s, vv = t, (62) 

whilst writing / = ln(u) for the potential corresponding to the choice of the parameters 
A = -1, B = C = 1, (58) becomes 

„ (s-a)(as-l) A (t-p)(pt-l) 

uu — T~\ ? uu — 77^ 7\ • (63) 

s(a 2 - 1) ' t(P 2 - 1) y J 
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System 


Parameters (A,B,C) 


Potential 


Equation 


Fi 


(-1,1,1) 


/ = ln(u) 


A 9* 

AZ p=2(a 2 + l)/(a 2 _i) 




(-1,1,-1) 


/ = 1b(«) 


A 9* 




(1,0,0) 


X 1 / \ 

/ = ln(u) 


A2 


Fn 


(0,1,1) 


/ = ln(w) 


H3* =4/a 2 




(0,0,1) 


/ = ln(w) 


H3 




(1,0,0) 


/ = u 




Fin 


fo 1 1) 


/ = « 


A lp=«2 




(1,0,0) 


/ = Hv) 


113,5=0 




(0,1,0) 


f = v 


A15=o,q^ q 2 


Fy 


(0,1,0) 


f = u 






(1,0,0) 


f = v 


HI 



Table 3. Equations governing different potentials of the systems (54)-(57). The 
potentials are introduced through equations (58)-(61) with the indicated choice of 
parameters (A,B,C). It is explained in the text how to construct a transformation 
between equations (idolons) that govern different potentials of the same system. 
Parameter associations in the listed equations are given connecting a and p or 
transforming a, similar associations connecting f3 and q or transforming f3 are implicit, 
they are omitted from the table for brevity. 



Eliminating s and t from (54) by means of (62) we find that v satisfies equation A2 (41), 
whilst eliminating s and t from (54) using (63) we find that variable u satisfies equation 
A2* (21) with the parameter associations 

a 2 — 1 p l — 1 

Equations governing the potential / for different choices of the parameters A, B and C 
are referred to in [5] as idolons of the underlying system governing s and t, in particular 
we have shown that equations A2 and A2* are idolons of (54). 

A Backlund transformation can be obtained by composing these relations, that is 
by eliminating s and t from formulas (62) and (63), 

(vv — a)(avv — 1) „ (vv — B)(Bvv — 1) 

uu = y - >\ \ uu = y - yjf ' . 64 

vv(a 2 — 1) vv(p z — 1) 

This is the Backlund transformation between A2 and A2* that appears in the first part 
of Table 2. The second transformation between these models can be obtained similarly 
after observing that the potential corresponding to choice of parameters A = C = — 1, 
B = 1 is also governed by model A2*. 

Table 3 contains the data required to construct, by the same procedure, all but the 
first two transformations in the first part of Table 2. 
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7. 5. Discriminant properties of the transformations 

All transformations listed in Table 2 are in the same class, they involve equations on 
lattice edges and the defining polynomial B = B(u, u, v, v) is degree-one in each of u and 
u, and degree two in each of v and v. This defining polynomial also has the following 
discriminant properties, 

{d u B)(duB) - {d u duB)B oc nfih(v,v), 

(d d B) 2 - 2{d\B)B oc r] 2 h*(u,u), 

where /j, and r\ are polynomials in v, and h, h* are the edge biquadratics associated 
with the two equations connected by the Backlund transformation. Specifically h 
coincides with Hi associated with the multi-affine equation in v through the generalised 
discriminant (46), while h* coincides with H 1 associated with the multi-quadratic 
equation in u through the discriminant formula (16). This discriminant property is an 
important feature of the transformations that combines asymmetrically the underlying 
discriminant characterisations of the multi-quadratic models (18)-(24) and the multi- 
affine ABS equations (38)-(45). In particular it may be used as a basis for their 
construction. 

One approach to such construction is offered by recognising the biquadratics 
themselves (6)-(13) can take natural discriminant forms, as the primary example we 
recognise that (7) is proportional to 

(u + u + 25a 2 ) 2 - A[a(u + 5)] [a(u + 5)}, 2a 2 = p + 1. 

It is easily seen that this expression emerges as the discriminant, with respect to v or 
v, of the expression 

[a(u + 5)}v 2 - [u + u + 25a 2 ]vv + [a(u + 5)}v 2 , 

which itself coincides with the polynomial defining the first transformation in Table 2. 
The precise choice of this expression is not unique, and in particular the biquadratic 
obtained from the generalised discriminant in u and u emerges a posteriori after choosing 
the expression on the basis that the required discriminant property in v and v is obtained. 

We remark that the classification of polynomials B with the discriminant properties 
(65) is interesting from the point of view of exhausting transformations in the same class 
as those listed in Table 2 (such task is similar to the one solved for multi-affine quad 
equations in [26]). 

8. Discussion 

The discriminant factorisation property, as laid out in Section 2, allows reformulation 
of the multi-quadratic quad equation as a system that defines single-valued evolution 
from initial data. Our main result is the list of models constructed on the basis 
of this hypothesis in Section 4, it contains all previously known integrable multi- 
quadratic quad equations as well as a substantial number of new models that exhibit 
the same integrability features. The nature of the relationship between the discriminant 
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factorisation property and the integrability is therefore an important question. Here 
we have made additional assumptions beyond the discriminant factorisation hypothesis, 
therefore although our investigations are suggestive, they leave open the problem of 
determining whether this property is sufficient for integrability in this class of models. 

Beyond elliptic (and hyperelliptic) functions, degree-two equations studied 
systematically from the point of view of integrability seem to be relatively rare. 
An isolated precedent exists within the framework of Painleve analysis for ordinary 
differential equations. A list of degree-two counterparts of the Painleve equations was 
obtained by Chazy [27], loosely speaking these are second-order ordinary differential 
equations with the Painleve property that are quadratic in the second derivative term. 
The difficult step of rigorously classifying this class of equations was made by Cosgrove 
and Scoufis [28]. In that setting one of the principal features is that the higher degree 
models do not define new transcendents (see [29]), in particular they are solvable in terms 
of the degree-one Painleve equations. Here the primary model, namely the equation that 
we identify as a multi-quadratic counterpart of Q4, is the only new model we have not 
yet connected back to an equation from the multi-affine class. The parallels between 
integrable quad-equations and the Painleve-type equations suggest that a Backlund- 
type transformation establishing such connection should exist. Verifying or falsifying 
this is therefore an important open problem. 
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